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A variety of Pioneer Venus observations suggest a global scale, day-to-night Venus ther- 
mospheric circulation. The two-dimensional hydrodynamic model of Dickinson and Ridley (1977) 
correctly predicted the gross characteristics of this largely symmetric circulation. However, it failed 
to calculate the observed cold nightside temperatures, and the exact phases and densities of the 
neutral constituents. This paper presents model studies of the dynamics and energetics of the 
Venus thermosphere, in order to address new driving, mixing and cooling mechanisms for an im- 
proved model simulation. 

The adopted approach has been to re-examine the circulation by first using the previous two- 
dimensional code to quantify those physical processes which can be inferred from the Pioneer 
Venus observations. Specifically, the model was used to perform sensitivity studies to determine 
the degree to which eddy cooling, eddy or wave drag, eddy diffusion and 15 gm radiational cool- 
ing are necessary to bring the model temperature and composition fields into agreement with 
observations. Three EUV heating cases were isolated for study. 

Global temperature and composition fields in good agreement with Pioneer data were obtained. 
Large scale horizontal winds < 220 m/s were found to be consistent with the observed cold 
nightside temperatures and dayside bulges of O, CO and C0 2 . Fine tuning required that an eddy 
coefficient <20% of previous one-dimensional models be used for nightside diffusion (K = 7.5xl0 6 
cm 2 /s). Very little eddy diffusion was required for the dayside (K<4xl0 6 cm 2 /s). Observed dayside 
temperatures were obtained by using a 7-19% EUV heating efficiency profile. The enhanced 15 j«m 
cooling needed for thermal balance is obtained using the best rate coefficient (K cc , 2 _ 0 = 5xlO -13 
cm 3 /s) available for atomic O collisional excitation of CO 2 (0, 1, 0). Eddy conduction was not 
found to be a viable cooling mechanism due to the weakened global circulation. The strong 15 g m 
damping and low EUV efficiency imply a very weak dependence of the general circulation to solar 
cycle variability. Finally, the NCAR terrestrial thermospheric general circulation model (TGCM) 
was adapted for Venus inputs using the above two-dimensional model parameters, to give a three- 
dimensional benchmark for future Venus modelling work. 
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Figure 3. Global mean Venus one-dimensional thermospheric temperature 
profiles: eddy cooling tests. Here K t is the turbopause value of the eddy 
coefficient and Z t corresponds to the turbopause level. 
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Figure 4. Venus two-dimensional model temperature response to eddy cooling. HSYM corresponds to the 
Keating et al. (1984) model atmosphere, and VTS3 to the Hedin et al. (1983) empirical model. 
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Figure 5. Venus two-dimensional model response to viscous drag tests. The effects 
are shown for various adjustments of a Rayleigh friction coefficient (K RA ) (a) 

Too, (b) terminator u(z) profiles. 
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Figure 6. Venus two-dimensional model 
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Figure 7. Family of terminator horizontal wind profiles. Various drag profile shapes and K RA 
magnitudes were examined to obtain the best day-night distribution of model composition and 
temperature fields. 
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Figure 8. One-dimensional radiative transfer model reference 
temperatures: dayside and nightside mean profiles. These were 
calculated using the enhanced 15 /tm cooling (K C02 _ 0 = 5xlCH 3 cm 3 /s) 
possible via atomic Ocollisional excitation of C0 2 (0, 1, 0). Results are 
given for the BEST EUV case (7-19% profile) heating. 
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Best Parameter Adjustments for Each EUV Efficiency Case 
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Figure 9. Best parameter adjustments for each EUV efficiency case. The 7-19% case was isolated for further study since the 15/mi 
cooling K C q 2 _ 0 rate adopted is comparable to that inferred by Sharma and Nadile (1981) from a terrestrial rocket measurement 
of C0 2 emission limb radiance above 100 km. 




Figure 10. Venus two-dimensional model temperature: BEST EUV case. Comparisons are made with 
Pioneer Venus model atmospheres (HSYM and VTS3) as well as with the adopted nightside 
temperature profile of Massie et al. (1983). Solar minimum fluxes of Heroux and Hinteregger (1978) 
(F10.7 = 80xl0~ 22 Wnr* 2 Hz’ 1 ) are used to derive a new noon temperature profile. Exospheric 
temperatures are seen to vary by nearly 70 °K from solar minimum to solar maximum conditions. 
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Figure 11. Heating and cooling rates at 60° and 120° SZA: BEST EUV case. Weakened adiabatic 
heating and horizontal advection on the nightside (plus improved 15 fim cooling) combine to give 
cooler nightside temperatures than Dickinson and Ridley (1977). 





Model and Empirical Data Set Comparisons at 150 km 
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Figure 12. Model and empirical data set comparisons at 150 km. 
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Figure 14. Solar minimum case contour plots (SZA vs Z): 

(a) Total temperature in deg K. 

(b) Horizontal velocity in cm/sec. 
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Figure 15. 60° SZA model helium (He) calculations. It is apparent from comparison of the one-dimensional 
and hydrodynamic model predictions for He at 100 km (Z ~ -8) that the self-consistent large-scale dynamics is 
crucial to a proper calculation of composition. One-dimensional models are insufficient. 




io " 5 icr 4 to -3 io ' 2 10 1 io° 

FRACTIONAL NUMBER DENSITY 


Figure 16. Venus two-dimensional model composition, helium: BEST EUV case. Noon and 
“nightside peak” profiles of helium are shown. The two-dimensional symmetric model ap- 
proximately simulates dynamical conditions for the observed 3-4 AM helium peak near LT 
= 23. 
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Figure 17. T + (u, v) contours on a pressure level slice for the Venus 
three-dimensional model adaptation of the terrestrial TGCM (a) Z = 
5, (b) Z =0, (c) Z = -4. Temperature contours are superimposed 
upon horizontal wind vectors. 
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Figure 18. Equatorial slices (SZA vs Z) for the Venus three 
dimensional adaptation of the terrestrial TGCM: 

(a) T (total temperatures) 

(b) h (height field in km) 
















